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§ ! Time-resolved optical lineshapes are calculated using a second-order inhomogeneous cumu- 

lant expansion. The calculation shows that in the inhomogeneous limit the optical spectra are 
determined solely by two-time correlation functions. Therefore, measurements of the Stokes- 
shift correlation function and the inhomogeneous linewidth cannot provide information about 
the heterogeneity lifetime for systems exhibiting dynamic heterogeneities. The theoretical re- 
sults are illustrated using a stochastic model for the optical transition frequencies. The model 
rests on the assumption that the transition frequencies are coupled to the environmental re- 
laxation of the system. The latter is chosen according to a free-energy landscape model for 
dynamically heterogeneous dynamics. The model calculations show that the available exper- 
imental data are fully compatible with a heterogeneity lifetime on the order of the primary 
^ ' relaxation time. 
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I. Introduction 



c3 

a 

^ ■ The non-exponential response of supercooled liquids and other systems exhibiting slow re- 
laxation has been a topic of many investigations, cf. [|TJ, 0. At present, there seems to be 
consencus about the fact that the primary or a-relaxation in supercooled liquids near the 
calorimetric glass transition temperature T g is of a heterogeneous nature^, £|. Several exper- 
imental techniques have been invented in order to monitor various aspects of these dynamic 
heterogeneities ||. A supercooled liquid above T g represents an ergodic system on long time 
scales. Therefore, the lifetime of the dynamic heterogeneities is finite. Up to now a cou- 
ple of experimental methods exist that are able to detect these lifetimes. All of them are 
based on the idea of selecting a sub-ensemble of slowly relaxing entities and monitor its 
re-equilibration afterwards. The development of such techniques started with the invention 
of reduced four-dimensional (4d) NMR|J, which subsequently has been applied to a number 
of glass-forming systems 0, || [| |10|| . Other relevant techniques are the deep-bleach experi- 
ment designed by Cicerone and Ediger[ll] and the atomic force microscopy method used by 



Russell and Israeloff JT2| . 

A different experimental technique to monitor the dynamical behavior of liquids is pro- 
vided by time-resolved fluorescence and phosphorescence spectroscopy, for a recent review 
see[13|. After an initial optical excitation with a laser-pulse, the time- dependent emission 
spectra are observed, cf. Fig.l. Most theoretical investigations of time-resolved fluores- 
cence have concentrated on the Stokes-shift correlation function, which is related to the 
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time-dependent mean (uj(t)) of the emission spectrum: 



C(t) 



(u(t)) - Hop)) 
(w(0)> - (w(oo)> 



In particular, the dependence of C(t) on molecular variables of the solvent have been inves- 
tigated for both, polar |I3|, |Tj| [H], and nonpolar solvation ]T9|, [2"0fl , see also the list of 
references in ref . ]13| . 

In the field of supercooled liquids, triplet state solvation experiments have been applied. 
The observed behavior of the time-resolved inhomogeneous broadening of the optical lines has 
been interpreted in terms of dynamic heterogeneities pi], |22|, |23|. For one well studied system, 
quinoxaline (QX) in supercooled 2-methyltetrahydrofuran (MTHF), it has been concluded 
that the lifetime of the dynamic heterogeneities or the 're-equilibration time', (r req ), exceeds 
the mean time scale of the a-relaxation, (r a ), considerably 0, |25 ]. This is a somewhat 
surprising result for the following reasons. In all 4d-NMR experiments, performed on a 
variety of glassformers, ranging from polymers0 to simple molecular liquids ||, (r req ) ~ (r a ) 
has been found. Furthermore, a strong temperature dependence of (r req ) has been observed 
in the deep-bleach experiments. These techniques cover the temperature range in which 
the solvation experiments have been performed. Therefore, MTHF is assumed to behave 
different from all systems studied so far. In addition, the finding of a long {r req ) relies on a 



purely phenomenological interpretation of the data in terms of a very special model[23. 25]. 

It is the purpose of the present paper to present a sound theoretical description of time- 
dependent optical line shapes in complex systems. The actual calculation is based on a so- 
called inhomogeneous cumulant expansion [P^]. The results obtained bear a formal similarity 



to those of ref.pBfl, thus partly validating the phenomological approach used there. In order 
to illustrate the general results, I then treat the optical transition frequencies as a non- 
Markovian stochastic process and calculate the relevant quantities in terms of a free-energy 
landscape model for the primary relaxation [27|, |29], [3(J. The outline of the paper is 



as follows. In the next section I will give a theoretical description of optical lineshapes in 
the inhomogeneous limit in terms of an inhomogeneous cumulant expansion. Some formal 
aspects of the derivation are treated in Appendix A. Sect. Ill gives a brief description of 
the free-energy landscape model which then is used to calculate the relevant quantities. The 
results are discussed in Sect. IV and the paper closes with some conclusions. 



II. Theory of time-dependent optical lineshapes 

In this section the theory of time- dependent optical lineshapes in the inhomogeneous cumu- 
lant expansion will be presented. Throughout the discussion the inhomogeneous limit will 
be considered exclusively and vibronic regressions and also the finite lifetime of the excited 
state will be neglected. This means, I focus on the time evolution of the (slow, (r a ) longer 
than roughly 10 _8 s) solvent relaxation. To this end, the following simplified Hamiltonian of 
a two-state electronic system is considered: 

H = \g)[H + W g ](g\ + \e)[H + W e + hu eg ](e\ (2) 

Here, the electronic ground and excited state are denoted by \g) and |e), respectively. H is 
the solvent Hamiltonian in the absence of the solute, W g and W e denote the solute-solvent 
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interactions in the respective electronic states of the solute and hu eg is the — transition 
frequency of the solvated solute. In the following, % will be set to unity. In Appendix A 
the time-dependent optical lineshape is calculated from the third order (in the electric field) 
response |17], [18] in the rotating wave approximation. In the inhomogeneous limit both, time- 
dependent fluorescence (|e) = \Si)) and phosphorescence (|e) = |T\)) can be treated on the 
same footing. In the latter case it has to be assumed additionally that intersystem crossing 
and internal conversion are much faster than any other relevant time scale. For systems 
with slow environmental relaxation with typical relaxation times in the ms ■ ■ ■ s regime this 
provides no severe restriction. 

In this section the optical emission lineshape will be calculated within the framework of 
an inhomogeneous cumulant expansion [^] starting from 



I(u,t) = du P(uj } t\uJo)Sg(uj ) (3) 

J — oo 

Here, S 9 (lj) denotes the initial absorption spectrum, usually approximated by the initial 
(non-equilibrium) emission spectrum. Additionally, P(u},t\u}o) is the conditional probability 
for finding the transition frequency u at time t given it had the value ujq at time t = 0. 
P{uj,t\ujo) can be cast in the form 



OO POO 
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P(uj,t\uj ) = - — , r dn dT 2 e- luJ0T1 e iuJT2 (e iuj(0)T1 e- tuj{t)T2 ) e (4) 

47T 2 Se(^ ) 



which follows immediately from the definition of the joint probability, P(u,t\LJo)S e (LjQ) = 
{S(uj — u(t))5(u! — u(0))) e with (A) e = Tr[Ap e ], cf. Appendix A and ref.pOQ. In Appendix 



A, Eq.([|) is derived from the general expression for the third order response and the in- 
homogeneous cumulant expansion is performed in detail, yielding Eqns. flAT5| , |A.7| ) for the 
lineshape. 

Eq. (|4]) allows to discuss the formal connection to observables of other experimental tech- 
niques such as 2d-NMR. In this case, the conditional probability P(ujl, t\uj' L ) with ojl denot- 
ing the Larmor frequency is observed direct ly|5I]. In contrast to the optical case in NMR 



the dependence of ujl on molecular orientation is known explicitly and equilibrium correla- 
tion functions of the form (e luJL ^ Tl e~ lu}L ^ T2 ) are accessible directly. As is well known |31 



2d-NMR spectra are not able to provide information about dynamic heterogeneities let alone 
their lifetimes. This is because in order to address this question one has to probe the system 
at four points in time at least |H| . It was this knowledge that led to the invention of higher- 
dimensional NMR techniques 0. It is important to note that according to Eqns. (|3],|3]) I(co,t) 
does not contain any information about molecular motions beyond the two-time correla- 
tion functions (e 1 ^^ 1 e~ m ^' T2 ) e and therefore cannot provide further insight into the issue 
concerning the lifetime of dynamic heterogeneities. Furthermore, due to the complex depen- 
dence of u(t) on all relative distances and orientations of the solute relative to the solvent 
molecules, the optical lineshapes can be interpreted in terms of molecular quantities in a sim- 
ple way only in favourable cases. An example is the relation between the solvation coordinate 



and the dielectric properties of the solvent in case of polar solvation |T3|, [L4|, [15], [16], |T7 |. 

Now, I proceed to perform an inhomogeneous cumulant expansion of Eq.(Q). For this 
purpose the Taylor expansion of the phase correlation function ( e M0)n e -M*) r 2^ e j s congjcl. 
ered. Setting (u) e = and a\ = (uj 2 ) e without loss of generality, one has in a first order 
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approximation 

(e M0)n e — (t)r 2)e „ x _ ( (T 2 / 2)( r 2 + r 2) + a e 2 C(t)nr 2 

with the Stokes-shift correlation function given by: 

C(t) = (co(tM0)) e /a 2 e (5) 

The non-exponential decay of C(t) can conveniently be parametrized via a Kohlrausch func- 
tion, C(t) = exp [(— {t/T K Y K )}. This can be rewritten as 

C(t) = J d\p{X)e- xt (6) 

with relaxation rates A and a distribution function p(X). One then has, using JdXp(X) = 1, 
( e Mo)n e -M*te) e ^ /dAp(A){l-(<7^/2)(rf+T|)+^e -A *rir2}. Exponentiating the expression 
in the curly brackets gives: 

( e «*(o)n e -Mt)^ e „ | rfAp(A) exp |_^ [r 2 + r 2 _ 2e - A V 1 r 2 ]| (7) 

It is illustrative to compare this result with the corresponding one of an ordinary cumulant 
expansion. In this case one finds, without using Eq.(||) 

/ e iw(0)-n e -iw(t)7a\ ~ e -(^/2)[(r 1 2 +r 2 2 -2r 1 r 2 C(i)] 

Therefore, this expansion does not take into account the inhomogeneities reflected in Eq.@. 
In the same way as an expansion of the form C(t) ~ / dXp(X)(l — Xt) ~exp [— (/ dXp(X)X)t] only 
yields an exponential decay the ordinary cumulant expansion can only produce Gaussian 
probability distributions. The inhomogeneous cumulant expansion instead allows one to 
account for a more complex behavior ||26||. 

Performing the Fourier transform in Eq.@, using Eq.(|7|), yields for the conditional prob- 
ability: 

P(u,t\wo) = J dXp{X)P w (u,t\uJo) 

1 / [u - uj e~ xm 




Note that Eq.(P) is the same as Eq. ( |A.8| ) for o- 2 {Y) = u 2 independent of the environmental 
variable T, which in Eq.(|^) corresponds to A. It is important to notice that both, the distri- 
bution p(A) and the decay rates A are determined via the Stokes correlation function C(t). 
In particular, the decomposition of P(u;,£|u;o) into the functions t\uo) is meaningful 

only if the A are static quantities. This means that the influence of possible exchange pro- 
cesses has to be incorporated in the definition of these rates. Thus, in exchange models that 
start from bare relaxation rates which are modified due to exchange, the A are the effective 
relaxation rates and p(X) the effective distribution. This will be further discussed below. 
Finally, the time independent absorption and emission spectra are given by: 



SJu) = — r = e -( w - A )V(a^) an d SJu) = p eq (u) = -^=e~ u /(2a ^ (9) 
/2^| F J^ 2 
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where A = (ui) g denoted the overall red-shift. Additionally, a 2 = ((u — A) 2 ) g , a 2 = (u) 2 ) e 
and (uj) e = has been used without loss of generality. In Eq.(|9]) the spectrum S e (u) is 
determined by ensemble averages with respect to the excited-state Hamiltonian, because the 
system occupies this state before photon-emission. 

Using the result of the inhomogeneous cumulant expansion, Eq.(||), in Eq.(|3]) gives a 
superposition of Gaussians for the optical lineshape I(u,t). The expectation value of any 
dynamical variable A can then be calculated according to 

(A(t))= JduA(u)I(u,t) = J did JdLU A(Lu)P(uj,t\LU )S g (LU ) (10) 

From this one finds the n th frequency moment using A = uo n . The time-dependent mean is 
directly related to the Stokes-shift correlation function, 



(u(t)) = A x C(t) = A x / d\p(\)e~ xt (11) 

and the second moment to C{2t): 

(oo 2 (t)} = a 2 + (A 2 + a 2 - a 2 e )C(2t) where C(2t) = J d\p(X) e - 2Xt (12) 

Therefore, the time-dependent variance is given by: 

° 2 (t) = a 2 e + (a 2 - a 2 )C(t) 2 + A 2 [C(2t) - C(t) 2 ] (13) 

This expression was found to give excellent agreement with the experimentally observed 
linewidth in the triplet state solvation studies on supercooled liquids^, ^l], [24H . 

The results for the optical lineshape and the frequency moments obtained above and in 
Appendix A can be reduced to expressions given earlier in the literature. If there is only one 
value for the environmental variable V (or A) Eq. ([A.5|) reduces to the expression already given 



above|2(| |33|. Note that C(t) does not necessarily decay exponentially in this case. This 
result is obtained by performing an ordinary cumulant expansion, cf. the discussion above. 
In that case the variance is given by <J 2 (t) = a 2 + (a 2 — a 2 )C(t) 2 and changes monotonously 
from (x 2 (0) = a 2 to cr 2 (oo) = a\ because of C(0) = 1 and C(oo) = 0. 

The result for the Pf\}(u, t\uo) derived above, Eq.(|J), can be mapped onto an expression 



given by Richert[f2l| if his loosely defined 'local response functions' r) are replaced by 
e~ xt . This point will be discussed further below. Of course, in any theoretical calculation 
of a conditional probability one can only obtain approximations to the two-time correlation 
function ( e M ) r i e -M*)T 2 ^ j n a se cond order cumulant expansion the information is further 
reduced to the autocorrelation function (u(t)u(0)). As mentioned already above, no two- 
time correlation function is able to provide any information about the lifetime of dynamic 
heterogeneity. 



III. A free-energy landscape model for the optical lineshape 

In order to illustrate the results derived above in this section the transition frequencies u(t) 
are treated as a non-Markovian stochastic process. For this purpose I consider the description 
of u(t) in terms of a free-energy landscape model. This approach has been used previously to 
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describe several aspects of the a-relaxation in supercooled liquids^, p5[ . The basic idea is 



the following. The a-relaxation is assumed to be associated with activated transitions among 
an extensive number of metastable glassy states or valleys in the free-energy landscape ]54j]. 
This activated dynamics is described as a stochastic process e(t). For simplicity, e(t) is 
modeled as a stationary Markov process. The corresponding master equation (ME) for the 
conditional probability (Green's function) G(e, t\e Q ) = Pi|i(e;£|e ) then reads f35fl: 

— G{e,t\e ) = Jde'K(e\e')G(e',t\e ) - K{e)G{e,t\e ) with K {e) = Jde'n(e'\e) (14) 



Here, ft(e|e') denotes the rate for a e' — > e-transition. The dynamics described by Eq.([14|) 
is supposed to be relevant for the environmental relaxation of the system. Of course, such 
a model can be applied to any system exhibiting activated dynamics, e.g. proteins, if the 
«(e|e') are chosen appropriately. In the following calculations I choose a simple globally 
connected model[^7| for the kinetics. After an escape out of valley (initial state e') with an 



activation energy (Ea — e') any other valley (final state e) can be reached with a probability 
given by its density of states, 77(e) oce~ e / 2cr(e ) : 

K{e\e') = kn-n&ef* (15) 

Here, is a prefactor and (3 = [kBoitzmannX^ 1 ■ In ref.|27j other models have been considered 
also. 

The main assumption of the model is that if x(t) is a process observed experimentally, 
then x(t) is allowed to change its value solely due to the e' — > e transitions. Therefore, x(t) 
can be considered as a 'slave process' of e(t)j36|. In particular, x(t) is not Markovian. In this 
context it is to be mentioned that in terms of Heuers Q parameter [Q, which is a measure 
of how long a relaxation rate remains correlated to earlier values, this model corresponds to 
Q = l by definition, i.e. {r req ) = (r a ). Previously, mainly the orientation of tagged particles, 
x(t) = f2(i), and the position, x{t) = r(t), have been considered |27|, p8fl . In the present 
paper, x(t) is identified with the optical transition frequencies, x(t) = uj(t). In order to 
describe the dynamics of uj(t), different approaches can be used. 

A. Langevin equations 

A simple approach is to assign a value w t (t) to each of the valleys e and consider the coupled 
Langevin equations: 

Cj e {t) = -K(e)uj e (t) + J de'K(e\e')LO £l (t) + L e {t) (16) 

where L e (t) is a delta-correlated white noise, (L e (t)) =0, (L t {t)L e i{t')) =Y e ^5{t — t'). For a 
finite number of valleys, the integral in Eq.(^) is to be replaced by a sum. After diagonal- 
ization of the matrix defined in Eq. ([16]) one has the decoupled Langevin equations 

dj x (t) = -Xuo x (t) + L\(t) (17) 

where A are the eigenvalues. These equations are directly related to Gaussian conditional 



probabilities |35|, P7|h determined by the correlation functions (u\(t)ui\(0)) = cr e e . There 



fore, these Langevin equations yield Eq.(RI) without further approximations. Due to the 
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structure of Eg. (|16|) the p(X) are determined by the G(e, t|eo) , cf. Eg. (|Lj|) . In an actual 
calculation, one needs to specify how the transition freguencies u e (t) depend on the envi- 
ronmental variables e. A very simple choice would be to neglect any correlation between 
cu(t) and e(t)|38[]. Without going into details of an actual calculation here, it suffices to note 
that for any such 'Langevin model' the inhomogeneous cumulant expansion is exact. This 
is because in a Langevin approach one starts from Gaussian fluctuations from the outset 
and u(t) represents a multivariate Ornstein-Uhlenbeck process|J^]. This means that u(t) 
is a non-Markovian Gaussian process. As a conseguence, Eg.(||) results for the conditional 
probability. 



B. Master equations 

Another approach allowing for a richer dynamics is to consider the ME for the composite 
Markov process {u(t),e(t)} as has been done previously p7|, |28fl . As will be shown below, 
in contrast to the Langevin eguations considered above, in such an approach uj(t) is non- 
Gaussian in general. In the ME, the transition rates k(uj, e\u', e') are chosen as simple 
products, k(uj, e\u', e') = K(e\e')A(uj\uj'). Here, the dimensionless functions A(u\u') determine 
the magnitude of changes in uj associated with a e' —>■ e-transition. Allowing only for very 
small freguency changes, u — >oj' = (j+5, 5<Cl, one can use a Kramers-Moyal expansion of 
the ME as described in Appendix B. As a result, the following eguation for the conditional 
probability is obtained: 

Pm(u,e,t\u ,e ) = -K(e)Pm(uj, e; t\u , e ) 



+ 



P A FP (uj)] J Gfc'/<e| e ')Pi|i(cu, e'; t\u , e ) 



where the Fokker Planck (FP)-operator A (uj) is defined by 



and 



"=2^ (20) 

The 'fluctuation parameter' p <ti 1 measures the scale of the cu-fluctuations relative to the 
width of the eguilibrium distribution. An analytic solution of Eg. (|l8|) is not feasible in 
general. However, one can perform an expansion in terms of the eigenfunctions of A fp (uj) 
which are Hermite polynomials if w (z)p7|: 

P Ml (u>, e; t\u , e ) = g H ^ H ^ Gn (e, t\e ) where = u"/2al (21) 



n =o 2 n n\ x /2na 2 



The Green's functions G n (e,t\e ) are the solutions of the rate eguations: 

G n (e, t\e ) = -«(e)G n (e, t|e ) + [1 - np\ J de'«(e|e')G n (e , , t\e Q ) (22) 

Note that Go(e, t\eo) = G(e, t\e ) according to Eg.(|l^). A comparison of these rate eguations 
with those occurring in the treatment of molecular reorientations shows that the term np 
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corresponds to Pi(cos8), the Z th order Legendre polynomial of the jump angle 9. A full 
solution of the problem requires the calculation of all G n (e, t|eo), n — 0, 1, • • • , oo. From 
this solution the stochastic process u{t) can be defined as a projection of the composite 
Markov process {u(t),e(t)}. For this purpose one has to integrate the joint probabilities 
over the variables e. In the special case of uncorrelated equilibrium probabilites, p eq (uj, e) = 



p eq (u)p eq (e), this yields the marginal distributions^? 



P(u,t\u ) = Jde p eq (e ) ydePi|i(u;,e;t|u;o,eo) (23) 
From this expression and Eqns.(^Ij |22|) it becomes clear that P(u,t\uo) can not be written 



as a superposition of Gaussians. This shows that even though Eq. (|l8l) resembles a Fokker- 
Planck equation, u{t) is not only non-Markovian but also non-Gaussian in general. 

Inhomogeneous cumulant expansion 

The non-Gaussian nature of the stochastic process u(t) in the approach starting from a 
ME allows to discuss the inhomogeneous cumulant expansion in detail. For this purpose 
I proceed in exactly the same manner as in Sect. II. There, it became clear that the only 
relevant quantity determining the conditional probability is the frequency autocorrelation 
function. Using the properties of the Hermite polynomials, this function can be calculated 
from Eq.(|2l|): 



(u(t)u(0)) = Jduj p eq (uj )uj JduuP(uj,t\uj ) = a\ Jde p eq (e ) JdeG^e, t\e ) (24) 
= al fde p eq (e ) fde /dA«$ (1) (e, A«)%)(e , A^e"^' 



Here, 3>(i)(e, A^) and ^(i)(eo, A^) are the right and left eigenvectors of the matrix defined 
in Eq.(p2|) for n = 1. The A 1 - 1 - 1 are the corresponding eigenvalues. Performing the same 
calculation as in Sect. II, one recovers Eq.(||) for the conditional probability, P(u,t\u ), with 

p(A«) = fde p eq (e ) /de$ (l3 (e,A«)* (1) (e ,A«) (25) 



Thus, the inhomogeneous cumulant expansion provides an explicit expression for the 'distri- 
bution of relaxation rates', p(A). 

Note, that a Gaussian approximation is equivalent to approximating the Green's func- 
tions G n {e,t\e ), n> 1, by G n (e,t|e ) ~ /dAi$ ( i)(e, A (1) )^ ( i)(e , A (1) )e-" A<1) *. Inserting this 



expression into Eq. (|2~1~D allows one to compute the sum over n resulting in a Gaussian ]37| . 
Using Eq.(]23|) yields Eq.(|]) with p(\) given by Eq. (|25|) . It is important to point out that 
the frequency autocorrelation function, Eq . (|24]) , also determines the conditional probability 
in an ordinary cumulant expansion, cf. the discussion in Sect. II. Therefore, in the present 
context Eq.(|^) is the result of an inhomogeneous Gaussian approximation. 

In general the assumption of Gaussian fluctuations presents an approximation, only. 
There are, however, two limits in which uj{t) becomes a Gaussian stochastic process. As 
shown in Appendix C a perturbation expansion of Eq. (|22|) yields a Gaussian probability 
distribution in the limits of short and long times, respectively. 

In the short time limit, one has to replace the rates A in Eq.(||]) by pn(e) and to set p(X) = 
p eq (e). The relaxation rates are therefore determined by the escape rates ft(e). This means 
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that any single e — > e' transition is effective for the relaxation of the transition frequencies 
u>(t). This limit therefore corresponds to a situation in which the u(t) fluctuate in a quasi- 
static environment. 

In the long time limit, only the average relaxation rate, p(n) = JdeK(e), determines the 
temporal behavior of P(u>, t\uj ) and one has to set p(X) — 8( X— p(k)) in Eq.(§). Therefore, in 
this limit u(t) is an Ornstein-Uhlenbeck process [|35| which means that u(t) is Markovian and 
Gaussian in this limit. Only the time scale for relaxation is set by the mean environmental 
relaxation rate (k). 

The inhomogeneous Gaussian approximation discussed above smoothly interpolates be- 
tween these two limits, if p^C 1. Here, the difference between an inhomogeneous cumulant 
expansion and an ordinary one can be seen most clearly. In both cases the correct long time 
limit is approached, as is well known from the theory of random walks [55[]. The correct short 
time limit, however, is only reached if the inhomogeneous cumulant expansion is performed. 



Model calculations 

In Fig. 2 the Stokes-shift correlation function C(t), calculated according to Eq.fl24]), and 
the time-dependent part of the variance, Eq. (pf) , are shown for p = 0.1 and a g = a e . The 
parameters, and cr(e), cf. Eq.(|l5|), are chosen in such a way that C{t) is well parametrized 
by a Kohlrausch function, C(t) = exp [— (i/ T ic)^ K ]- For the curves with j3 K = 0.5, the choice 
was guided by the experimental results of Wendt and Richert on triplet state solvation in 
MTHF[^4j. As already mentioned in the last section, the shown behavior for C(2t) — C(t) 2 = 
[a 2 (t) — a 2 }/ A 2 is in excellent agreement with the experimental data. Comparing the different 
curves in Fig. 2, it is seen that for more stretched C(t) the range of positive C(2t) — C(t) 2 
is increased. Also the absolute value increases with decreasing (3k- Without showing the 
results of calculations for other choices of the fluctuation parameter p here I just mention 
that this has only a small effect on the overall behavior if p is not chosen too large. In 
this latter case, however, the Kramers-Moyal expansion is no longer meaningful and Eq.flT8|) 
ceases to be valid. 

As has been pointed out, the solution of Eq . fll8|) results in a stochastic process uj(t) that 
is not only non-Markovian, but also non-Gaussian. A detailed study of the deviations from 
Gaussian behavior of the exact solution of Eq. fllSD is difficult in general. However, one can 
compare the low-order moments calculated in the inhomogeneous Gaussian approximation, 
Eqns.(|8l,|25|), with those obtained using the exact Eqns. (|21|j22] ). In a way similar to the 
calculation of {u(t)u(0)) one finds for a e = a g 

(u(t)) (NG) = AxC(t) and (u 2 (t)) {NG) =a 2 + A 2 Jdep e i(e)Jde G 2 (e,t\e ) (26) 

The first moment of course coincides with the Gaussian approximation, Eq. flTTl) , {uj(t))(NG) = 
(u>(t)). For the second moment and the variance, cr 2 (t), this is not true. Here, Eq.flTSP has 
to be compared to 

a 2 (t) {NG) = a 2 e +A 2 at) where £(t) = J dep e "(e) j ' de G 2 (e,t\e ) - C(t) 2 (27) 

This quantity is plotted in Fig.3(a) for the same parameters as used in Fig. 2 for (3k = 0.5 
along with C(2t) — C(t) 2 for comparison. It is clearly seen that exceeds C(2t) — C(t) 2 



9 



in the whole range where both functions are nonzero. In order to qualitatively discuss 
higher order moments (u n (t))^c) it is sufficient to notice that these are determined by 
Jdep eq (e)JdeoG n (e,t\eo), as can be seen from Eq. (f21|) . Therefore, the functions 



C„(*) = / depute) J de G n (e, t\e ) - C(nt) (2* 



give a qualitative measure for the deviations of (u> n (t))(NG) from a Gaussian behavior. This 
is because in the Gaussian approximation the ( n (t) vanish identically, ( n (t) = 0. The ( n (t) 
are plotted in Fig. 3(b) for n — 2, 4, 6, 8. From this plot it is evident that the deviations from 
Gaussian behavior even increase with n. Thus, if one wants to go beyond the Gaussian 
approximation, it is not sufficient to consider only the first few moments. In particular, the 
second moment gives only a poor approximation to the inhomogeneous linewidth. Addition- 
ally, the Stokes-shift correlation function is no longer determined by the first moment but 
also higher moments become important. However, it should be kept in mind that if higher 
order terms are kept in the cumulant expansion of P(u,t\uo), these terms also have to be 
taken into account for the static spectra S e (co) for reasons of internal consistency, cf. the dis- 
cussion in Appendix A. As will be discussed further below, it is the identification of a 2 (t)(NG) 
with the inhomogeneous linewidth for one special exchange model that let Richert |25|] to the 
incorrect conclusion that models with (r req ) ~ (r a ) are incompatible with the experimental 



data obtained for QX in MTHF|g4| 



IV. Discussion 

In Sect. II and Appendix A an inhomogeneous cumulant expansion has been performed to 
calculate the time-dependent optical linewidths in complex systems like supercooled liquids 
or proteins. Because the expansion has been truncated at second order, the transition 
frequencies u(t) constitute a stochastic process that is Gaussian but non-Markovian in this 
approximation. In Sect. Ill the theoretical results have been illustrated by a specific stochastic 
model for u(t). The coupling of these frequencies to the environmental fluctuations, treated 
in terms of the free-energy landscape model, renders u(t) non-Markovian. The same is true 



if exchange models are considered instead [ 3£ 



In theoretical treatments of polar solvation JT^, [16], [17] a simple interpretation of the 
time-dependent Stokes shift C(t) and accordingly the solvation coordinate is given in terms 
of the dielectric properties of the solvent. If the solvation coordinate can be modeled 
as a Gaussian Markov process, also the interpretation in terms of solvation free energies 
is straightforward [15, 17]. This, however, usually leads to exponential relaxation. If in- 



stead the Stokes-shift correlation function decays non-exponentially, such a simple relation 
is no longer obvious as in this case the solvation coordinate does not constitute a Markov 
process. Rewriting the non-exponentially decaying Stokes-shift correlation function as in 
Eq.(^) still allows to relate C(t) to the dielectric response via an expression of the form 
C(t) ~ JdXp(X) exp (— t / Tdiei{\)) ■ Here, Tdiei(X) denote dielectric relaxation times, the val- 
ues of which depend on the model considered |13|. Thus, the treatment of u(t) as a projection 
from a higher-dimensional Markov process along with the inhomogeneous cumulant expan- 
sion might be helpful in the interpretation of the solvation coordinate in more complex 
situations. 
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The free-energy landscape model for the dynamics used in the last section to treat the 
spectral diffusion process is of a phenomenological nature and specific assumptions about 
the dynamical evolution of the 'order parameter' e(t) and the coupling between the different 
processes have to be made in order to obtain a tractable model. One of the main features of 
this model lies in the fact that there is no difference between 'relaxation rates' and 'exchange 
rates' as in typical exchange models|^0|- Additionally, as mentioned above, the choice of the 
n{e\e') according to Eq.fllS]) corresponds to a minimum rate memory, Q = l[^]. This means 
that the the lifetime of the dynamic heterogeneities is the same as the time scale of the 
a-relaxation. 

The treatment of the spectral diffusion process in terms of coupled Langevin equations, 
Eq.fllBl), immediately gives the result of the inhomogeneous cumulant expansion, Eq.(H), 
for the conditional probability. Consequently, the first moment is given by Eq.flTTD and 
the time-dependent variance by Eq. (|T3|) . The reason for this behavior is given by the fact 
that in such a Langevin equation model one starts from Gaussian fluctuations from the 
outset. Therefore, uj{t) represents a non-Markovian Gaussian process. According to ref.p4| 
this means that this model gives excellent agreement with the experimental data on QX in 
MTHF if the transition rates are chosen in such a way that the Stokes-shift correlation is 
well parameterized by a Kohlrausch function with (3k = 0.5. Therefore, the experimental 
data are fully compatible with Q = l. 

A richer dynamical behavior of u(t) can be studied by considering the more general case 
of a ME for the composite Markov process {u;(t), e(t)}. The Kramers-Moyal expansion of this 
ME has led to the introduction of the parameter p = (5 2 /2a 2 ) that measures the fluctuations 
in transition frequency relative to width of the equilibrium steady state emission spectrum. 
It has already been mentioned in the context of Eq.(p2|) that the fluctuation parameter p 
plays a role very similar to the jump angle in models for the reorient at ional motion. 2d-NMR 
experiments performed in the time domain allow to extract these jump angles, roughly on 
the order of 10° in supercooled liquids]^, P, f£[|. As it is the reorient at ional and translational 
motion of the solvent molecules that is responsible for the (^-fluctuations, one expects that 
p also takes on finite values. Note that such finite fluctuations cannot be treated in terms 
of Langevin equations. 

In general, the process uo(t) is neither Markovian nor Gaussian in the treatment using the 
ME. Only in two limiting situations the result of the inhomogeneous cumulant expansion, 
Eq.(||), is found for the conditional probability for small p. At short times the transi- 
tion frequencies fluctuate in a quasi-static environment and every e — > e' transition gives 
rise to relaxation. Therefore, in this limit the relevant rates are just the 'escape rates' 
= jde' tv(e'\e). Conversely, at long times the environmental relaxation appears to be 
fast compared to the changes in uj(t). Thus, only the average rate (k) is relevant for the 
time scale of the ^-fluctuations, rendering u(t) an Ornstein-Uhlenbeck process. Note that 
these two limits can be compared to the so-called short-time approximation and long-time 
approximation in stochastic lineshape theory [Q. Apart from these limits, u(t) is not a 
Gaussian process and the deviations from Gaussian behavior have been discussed in context 
of Fig.3. Furthermore, I have discussed an inhomogeneous Gaussian approximation for the 
conditional probability which allows to smoothly interpolate between the limits of short and 
long times, where it becomes exact. For this approximation, the result of the inhomogeneous 
cumulant expansion, Eq.(§), is recovered. In contrast, an ordinary cumulant expansion fails 
to reproduce the exact behavior at short times. 
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It has been mentioned already that the expression for the P(X)(u},t\uJo), Eq.(H), can be 
compared to a corresponding one given by Richert[Q, if his x(t, T ) are replaced by exp (— Xt). 
However, in ref . |23| these functions were denoted as conditional probabilities and an expres- 
sion similar to Eq.(U) was postulated to hold rather than derived from a theoretical treatment. 
Additionally, the meaning of the 'local response function' r) remains unclear. Unfortu- 
nately, no definition for these objects was given apart from the fact that they are related to 
the Stokes-shift correlation function via C(t) = (x(t, r)) which contradicts the term 'response 
function'. Thus, it appears that a meaningful interpretation of the x{t, t) can be given only 
if these functions are identified with exp (—Xt). Surprisingly, this apparently is not Richert's 
interpretation because he concludes that the time-dependent variance is determined by the 
'spatially distributed local solvent responses' r) in case of a static heterogeneity. Addi- 
tionally, in ref.||25|| he considers exchange models in which the 'time constants r of the local 
responses' fluctuate in time. He furthermore claims that the experimental results are only 
compatible with slow rate exchange. Of course, in an inhomogeneous cumulant expansion 
the rates A (or corresponding time constants r) have to be static quantities. As I have shown 
in the last section, stochastic models which do not distinguish at all between 'relaxation' and 
'exchange' give results compatible with the available experimental data[24|. The same holds 
for exchange modelsPDJ. The confusion seems to have its origin in the fact that the exchange 



models studied in ref.[^5j are not Gaussian. Therefore, a Gaussian approximation has to be 
applied to the models in order to be able to compare the results to an inhomogeneous cu- 
mulant expansion truncated at the second order. Otherwise, C(t) is not determined by the 
first moment and the linewidth is not given by o- 2 (t). As noted above, including higher order 
terms is non-trivial, as this should also be done for the static spectra. These difficulties can 
be avoided if one considers models which are Gaussian intrinsically, like Langevin equation 
models. 



IV. Conclusions 

In the present paper I have derived an expression for the time-dependent optical linewidths in 
the inhomogeneous limit for systems exhibiting slow relaxation. For this purpose, an inhomo- 
geneous cumulant expansion has been performed, utilizing the fact that a non-exponentially 
decaying Stokes-shift correlation function can be represented as a weighted sum of exponen- 
tially decaying functions. 

In a free-energy landscape model developed to model the primary relaxation in super- 
cooled liquids, the optical transition frequencies were treated as a stochastic process u(t), 
which is not Markovian in general. Two variants of this model have been considered, a 
Langevin approach and a treatment in terms of a master equation. In the former case the 
model is solved exactly by the inhomogeneous cumulant expansion as the fluctuations truly 
are Gaussian. Therefore, in this case u(t) is a non-Markovian Gaussian process. This does 
not hold for the more general ME approach. Here, the conditional probabilities are Gaus- 
sians only in the short time and the long times limits. I have additionally considered an 
inhomogeneous Gaussian approximation that interpolates between these limits. In general, 
however, in such more complex models u(t) is not only non-Markovian but additionally 
non-Gaussian. 

Due to the fact that the system is probed at two times only in the optical experiments 
under consideration, one cannot obtain any information about the lifetime of dynamic het- 
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erogeneities from these experiments alone [3^]. This is in contrast to recent claims where 
it was concluded on the existence of long lived heterogeneities within a special model sce- 
nario that was assumed to apply [p5[]. The model calculations performed in the present paper 
clearly show that the experimental data are also fully compatible with the assumption that 
the lifetime of the dynamic heterogeneities is on the same order as the a-relaxation time. 
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Appendix A: Nonlinear response theory of optical lineshapes 



In this Appendix the result for the time-dependent optical lineshapes will be derived from 
an inhomogeneous cumulant expansion of the corresponding third order optical response. 

The starting point is the Hamiltonian, Eq.(fJ), to which the radiation- matter interaction 
responsible for the electronic transitions is added. In the Condon approximation the latter 
is simply given by Hint = \g)(e\ + \ e )(g\- Treating the excitation pulse classically and 
the emitted radiation quantum mechanically, the following expression for the time-resolved 



lineshape is obtained in the rotating wave approximation, see eg. |0, ITS 



OO POO 

3 i dr 2 / dr ie 4{ ^* 1+ws * 3) 

— OO J — OO J — OO 



S(u>s,(^L,t) oc Re / dr 3 

J — OO 

E L {t -h-t 2 - t 3 )E* L {t -t 2 - t 3 )R(h,t 2 , t 3 ) 



(A.1) 



where Ei,(t) denotes the time-dependent amplitude of the incoming field and R(ti,t 2 ,t 3 ) 
a third order response function which in the inhomogeneous limit is given bypj 



is 



R{h,t 2 ,t 3 ) =Tr 



-iUt 3 e -iC e t 2 e iUti 



(A.2) 



In this expression, U = W e — W g and g g denotes the canonical ground state density operator, 
g g = exp (—j3H g )/Tr[exp (—f3H g )} with H g = (g\H\g) and (3 = (ksT)^ 1 . Furthermore, C e 
denotes the Liouville operator corresponding to if e ||26|]. 

Eq.( |A.l| ) describes the processes visualized in the double-sided Feynman diagram p6 



in 



Fig. 4. It is to be noted that a term corresponding to a resonant Raman process has been 
neglected completely in Eq. ( |A.1| ). 



Following ref. |T7j], in a next step it is assumed that the solvation coordinate U is the only 
relevant quantity determing the solvent relaxation. After a projection onto the operator 
U and a second order cumulant expansion, one gets a simple expression for the response 
function R(ti,t 2 ,t 3 ), which is determined solely by the Stokes-shift correlation function 
Tr {e lHet Ue~ tHet }U e (0)g e and some static quantities. Therefore, in the inhomogeneous 
limit the optical lineshape is completely determined by a two-time correlation function. 

In order to perform an inhomogeneous cumulant expansion ||26||, it is assumed that all 
quantities depend on a configurational variable T chosen from some distribution p(T). One 
then calculates R(ti, t 2 , t 3 ) = j dTp{T)R{ti, t 2 , t 3 ; T), cf.|^B|, with a result analogously to that 
obtained in ref. [pjj : 



R(h, t 2 , t 3 ) = J dTp(T) exp (im e (r)fi) exp (it 3 [m e (T) + (m g (T) - m e (T))C(t 2] T)} 



x 



x 



where 



x exp (-^/2[a e 2 (r) + (a*(T) - a 2 e (T))C\t 2 ; T)]) 
x exp (-^/2^ 2 (r) - ht^JT^ih; T)) 



m a (T) = (w(0; r))„ , a 2 a (T) = ((lu 2 (0; T)) a - m 2 a (T) 
C(t; T) = «(«(*; r)«(0; T)) e - m 2 e (T)) / a 2 e (T) 



(A.3) 



(A.4) 



Here, (A(t;T)) a — Tr[A(t;T)g a ], a&{e,g} and u(t;T) is the classical analogue of (cj eg (T) + 
U{t] T)). As in Eq.(^|) one has (A(t)) a = fdTp(T)(A(t; T)) a . 
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A simple expression for S(us,uJi,t) is obtained by assuming that t\ is very short and 
additionally neglecting the optical coherence during t 3 . In this case t 2 can be replaced 
with the observation time t, cf. Fig. 4. Formally, these approximations amount to replace 
E L {t-t x -t 2 -tz)E* L {t-t 2 -h) in Eq. ([A~lD by |^| 2 5(t 2 -0 5 (*i)- The remaining ^-integration 
can be performed and one finds for the (^-independent normalized lineshape: 



r( * hv 1 ( Ws - K r + (m g r - m e r c t ; r } 

J (ws, t) = / dTp(r) - exp " 



2 - 



with 



E 2 (t; r) = a e 2 (r) + (a 2 (r) - a 2 (r))C 2 (t; T) 



9 

Alternatively, Eq. ([A.5 ) can be cast in the form 



I(u,t) = JdTp(T) Jdw'P {v) {u,t\u')S g {u']Y) (A.7) 

where S a (uj;T) = [2ttct 2 (T)}- 1 / 2 exp [(_(^ _ m a (T)) 2 / (2a 2 a (T))}) . P (v) (u,t\u') is the condi- 
tional probability to find the emission frequency uj at time t, given that it had the value u' 
at time t = 0, 

P (r) (u,t \u ) = exp - L° 1 ' 2 / J F ^ A ' 8 

/27ro- 2 (r) i-c 2 t ; r) V 2a 2 (r)(i-c 2 (t;r)); 



From the properties of the conditional probability one has 

l(u,0) = JdTp(T)S g (uj;T) and oo) = J dTp(T)S e (uj; T) (A.9) 

One might argue that the derivation given here is rather formal, especially as the result 
again can be cast into a form that resembles ordinary linear response theory, Eq. ( |A.7|) . The 



reason for the formal treatment is twofold. From the inhomogeneous cumulant expansion for 
the response function, Eq. (|A.3|) , it becomes clear that the 'ingredients' in Eq. (|A.7| ) are not 
independent. This does not become clear from Eq.( |A.7p alone. If only this expression would 
be given, there would be no reason for restricting the cumulant expansions for S g (u';T) 
and the conditional probability to the same order. Also, a generalization of the expression 



I(uJ,t) = rf^, duoP(u>, t\u> )Sg(uo), cf. ref . P0| , is not straightforward in general. In particular 
it is not clear how to perform the inhomogeneous cumulant expansion. Only in the limit of 
T-independent static quantities, it suffices to consider P(u,t\uo). 

In the text, only the case of T- independent static quantities m a and a a is considered for 
simplicity. This means that one has S a {u) = S a (uj;T) independent of T. Additionally, this 
allows to rewrite Eq. ([A.7|) as t) = /f^ du) P(u, t\u} )S g (ui ). In the text I choose m e = 
and additionally rename m g = A denoting the overall red-shift. 

Appendix B: Kramers Moyal expansion 

In this appendix the derivation of Eq. (|T8|) is outlined. The starting point is the ME 



J^Pi|i(cj, e; t\u , eo) = J du'A(uj\uj') J de'/c(e|e / )Pi|i(w / , e'; t\u , e ) (B.l) 

duo' k(uj'\uo) I de' K{t'\t)Pi\i{u , e; t\oj , e ) 
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As noted in the text, the A(ui\u') are assumed to be finite for small jump- lengths only. This 
means, that u(t) is modelled as a so-called one-step process^, |43|] with a jump length 5 
(5 2 = 2a 2 p, cf. Eq. fl2~0|)). Furthermore, it is assumed that spectral diffusion in the limit of 
vanishing 5 can be viewed as a diffusion in a potential of the form U(uo) = (4:f3<j 2 )~ l uo 2 with 
/3 = (/c^T)™ 1 . (The factor 1/4 is introduced here for convenience only.) The force originating 
from this potential ensures that the time-dependent distribution approaches the equilibrium 
distribution p eq (u), Eq.(|]), in the long time limit. With the functions A(u) = —puo and 
B(uo) = 2pa 2 the A(u\ to') can be written in the from 



A(u\u') 



25 



+ 



Biuj) 
25 2 



5{u' -[u + 5]) + 



A[u) 
25 



+ 



B_[uj) 
25 2 



5(w'-[u-S\) (B.2) 



The occurence of the inverse powers of 5 has its origin in the scaling properties of the 
A(a>|u/)[f|3]]. In the limit of small 5 a Kramers- Moyal (KM) expansion (or equivalently in 
this case van Kampen's system size expansion |p5||) can be performed. With the definition of 
the jump moments a n (uj) = Jdu'(ou' — uj) n A(uj'\u)) |35| and the relations 



duj'A(uj\J)f(uj r ) = 



n=0 




a n (u)f(uj) and / dw'A(u\uj')f(u) = a (u)f(uj) 



holding for an arbitrary function f(u) one finds 



d . , . 
— Pi|i(o;,e;t|a;o,eoJ 



-a (w)K(e)Pi|i(a;, e; t\u , e ) 



+ 



a (a;)+A^ M H / dt>K{e\(?)Pm{u, e'; t\u , e ) (B.3) 



UJ 



E 



n=l 



-i)" (_d_y 



a n [uj) 



Using ao(w) = 1, = —pu and 02(00) = p2a\ 



with k KM [ 

in Eq.( p.3|) and truncating the expansion at the second order directly yields Eq . (|i~8|) in the 
text. The validity of the truncation is controlled by the smallness of the jump moments 
a 3 (ou) = —2p 2 a 2 e u and 0:4(0;) = ^o\p 2 . 



Appendix C: Perturbation theory for /?<Cl 

This Appendix gives a brief summary of the perturbation expansion of Eq. (|22]) . For small 
values of the fluctuation parameter p = 5 2 / (2a 2 ) a first-order expansion yields 

G n (e, t\e )~G(e, t\e )-np J*dr J de' J de"G(e, t-r\e')n(e'\e")G(e" , r|e ) (C.l) 

and thus, using JdeG(e,t\e ) = 1 

J deG n (e, t\to) — exp —np J dr J de f K(e)G(e , r|eo) (C.2) 

Here, exp (— x) ~ (1 — x) and the definition of n(e), cf. Eq.([L4]), has been used. Insert- 
ing Eq. (|C.2|) in Eq.(|23|) and using Eq . (|2l|) shows that the resulting sum over n can be 
evaluated [0. 
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In the short time limit, one can replace G(e', r|eo) by its initial value, G(e',r = 0|eo) = 
8(e' — 6o), yielding JdeG n (e, t\eo) — exp [— npK(eo)t\). Thus, one finds Eq.(||) with A = pfc(e) 
and p(X) = p eq (e). 

In the long time limit, on the other hand, one can write G(e',r = oo|eo) = p eq (e') in 
Eq. flC.2| ). This gives / deG n (e, t\e ) ~ exp [— np(n)t], where the mean relaxation rate is defined 
as (k) = fden(e). Again, one recovers Eq.(H) with p(X) —8(X— p(k)). 

Figure captions 

Fig.l : Sketch of the evolution of the optical emission spectrum. The environmental re- 
laxation gives rise to the Stokes shift of the observed time-resolved emission spectra. 
Here \g) and |e) denote the ground and excited electronic state, respectively. I(ui,t) 
is the time-dependent optical lineshape evolving from the initial absorption spectrum 
S g (u) at t = to the steady state emission spectrum S e (u) at long times. 

Fig. 2 : C{t) (upper panel) and C(2t)— C(t) 2 (lower panel) versus scaled time t/rx for 
p = 0.1. The parameters are chosen such that the C(t) are well fitted to C(t) = 
exp {-{t/T K Y K ). Dashed lines: (5 K = 0.65, = 2.3 x 10 3 s _1 , a(e) = 3.055; full lines: 
(3 K = 0.50, = 7.2 x 10 4 s-\ a(e) = 3.65; dotted lines: (3 K = 0.35, = 5.83 x lO 9 ^- 1 , 
er(e) = 5.03. Here, cr(e) is given in temperature units ((3 = 1). Additionally, the values 
of e are restricted to the interval e G {— a(e) 2 , a(e) 2 }, giving equilibrium probabilities 
p eq (e) oc i](e)e~ l3t that are truncated Gaussians. Note, that the prefactors include a 
factor exp (— (3Ea) and therefore are only related in an indirect way to the attempt 
frequencies, typically on the order of 10 15 s _1 . 

Fig. 3 : (a): £(t) as defined in Eq.(p7|) and C{2t) — C(t) 2 versus scaled time for the same 
parameters as in Fig.2 for (3k = 0.5. The difference between £(t) and C(2t) — C(t) 2 
shows the deviations from the Gaussian approximation, Eq.(^), for the conditional 
probability. 

(b): The quantities ( n (t), defined in Eq.(^), versus scaled time. In the Gaussian 
approximation the ( n (t) vanish identically, ( n {t) = 0- 

Fig. 4 : Double-sided Feynman diagram for the excitation and for the emission of an elec- 
tronic two-level system, cf. Eq.( [A.l|) . The solute-system initially is in thermal equi- 
librium in the ground state (^>cx \g)(g\). The interaction with the incoming radiation, 
the wavy lines, produces an optical coherence (|e)(g|), which then is converted into an 
excited state population (|e)(e|). After photon emission, which again is accomponied 
with a coherence, the system is again in a ground state population (\g)(g\). 
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